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ABSTRACT 

The observed parameters of the young superstar cluster M82-A1 and its as- 
sociated compact HII region are here shown to indicate a low heating efficiency 
or immediate loss, through radiative cooling, of a large fraction of the energy in- 
serted by stellar winds and supernovae during the early evolution of the cluster. 
This implies a bimodal hydrodynamic solution which leads to a reduced mass de- 
position rate into the ISM, with a much reduced outflow velocity. Furthermore, 
to match the observed parameters of the Hll region associated to M82-A1, the 
resultant star cluster wind is here shown to ought to be confined by a high pres- 
sure interstellar medium. The cluster wind parameters, as well as the location 
of the reverse shock, its cooling length and the radius of the standing outer HII 
region are derived analytically. All of these properties are then conffimed with 
a semi-analytical integration of the flow equations, which provides us also with 
the run of the hydrodynamic variables as a function of radius. The impact of the 
results is discussed and extended to other massive and young superstar clusters 
surrounded by a compact HII region. 

Subject headings: Galaxies: star clusters — galaxies individual (M82) — ISM: 
bubbles — ISM: HII regions — ISM: individual (M82-A1) 



1. Introduction 



In many starburst, in interacting and in merging galaxies, a substantial fraction of star 
formation is concentrated in a number of compact [Rsc ^ 10 pc), young (a few 10^ yr) 
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and massive ( 10^ < Msc < 10^) stellar clusters (superstar clusters; SSCs). These entities 
may represent the dominant mode of star formation in these galaxies (McCrady et al. 2003; 
Whitmore 2006) and if held by gravity they also represent the earliest stages of globular 
cluster evolution (Murioz-Tufion et al. 2004). Their large UV photon output and their 
powerful, high velocity gaseous outflows (the star cluster winds) are now believed to be the 
major responsible agents for the large-scale structuring of the interstellar medium (ISM) and 
for the dispersal of heavy elements within their host galaxies and the intergalactic medium 
(IGM). 

M82-A1 is one of these SSCs in the galaxy M82 (Smith et al. 2006), with a half-light 
radius Rsc = 3 ± 0.5 pc, an age, tsc = 6.4 ± 0.5 Myr, an UV photon output A^*°* = 
(7.5 ± 3.0) X 10™ s-^ and a mass, Mgc = (1-3 ± 0.5) x 10^ Mq if one assumes an IMF 
with a Salpeter slope and 0.1 Mq and 100 Mq as the lower and upper mass cutoffs. Stellar 
evolution synthesis models (such as Star burst 99; see Leitherer et al. 1999) predict an average 
mechanical energy input rate for M82-A1 equal to Lsc = 2.5 x 10^*^ erg s~^. Given such a 
large energy deposition rate, it seems surprising that M82-A1 is, after such an evolution time, 
still surrounded by a compact HII region, with a radius Rhii = 4.5pc, a density uhii ~ 1800 
cm~^, a metallicity ~ 1 — 2 times solar and a total mass of only 5000 Mq, instead of having 
produced a large-scale superbubble. Note that a similar argument can be made for many of 
the SSCs embedded in the 150 pc nuclear region of M82 (see de Grijs et al. 2001; Melo et 
al. 2005) as well as for a number of the extragalactic HII regions (Kewley & Dopita, 2002; 
Dopita et al. 2005). 

Smith et al. 2006 suggested that the HII region associated to M82-A1 may result from 
the photoionization of a 4.5 pc shell blown into the ISM by the central cluster. However as 
the observed size of the HII region is not consistent with the standard bubble model (see 
Weaver et al. 1977; Bisnovatyi-Kogan & Sihch 1995 and references therein), they assumed 
that the kinetic energy of the photoionized shell represents only a small fraction (~ 0.1) 
of the mechanical energy, Lsc, released inside the cluster and that the high pressure of 
the ambient medium has reduced the expansion velocity to 30 km s~^, the FWHM of the 
observed emission lines. 

This scenario, although intuitively correct, rises several questions. First, the spectro- 
photometric age of the cluster, tsc ~ 6.4 Myr (Smith et al. 2006), is inconsistent with 
the kinematic age of the shell (r^) even if one assumes 30 km s^^ as the expansion speed 
throughout its evolution, = Rnii/Vexp ~ 1.5 x 10^ yr. Second, if the matter injected by 
supernovae and stellar winds is accumulated inside the HII radius, Rhii = 4.5 pc, then the 
gas number density should be, uhii > 3LscTsc/'^'^I^V^R%jj ^ 2.7 x 10^ cm~^, where is 
the mean mass per ion, and a wind terminal speed = 1000 km s~^ was assumed. This 
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density is an order of magnitude larger than that found in the HII region. Note also that 
only a small fraction, fnu < 0.001, of this gas could be photoionized by M82-A1 and thus 
in this interpretation, a good fraction of the mass supplied by supernovae explosions and 
stellar winds is hidden somewhere within the 4.5 pc volume. These arguments imply that 
the cluster and HII region parameters are not consistent with the wind driven bubble model. 

In order to have a consistent model, here we look at the detailed hydrodynamics of the 
gas reinserted by winds and supernovae within the cluster volume. The original adiabatic star 
cluster wind model, proposed by Chevalier & Clegg (1985), described the various assumptions 
that have been used by all followers. Assumptions such as the equal spacing between sources 
within the SSC volume and the fact that random collisions of the ejecta from nearby stellar 
winds and supernova explosions, were to result into a full thermalization of the matter 
reinserted by the evolving sources. Under such conditions, the pressure of the thermalized 
matter would exceed that of the surrounding ISM and the injected gas would be accelerated 
to leave the cluster with a high velocity, while establishing a stationary wind. For this to 
happen the gas has to acquire a particular velocity distribution, increasing almost linearly 
with radius from the stagnation point, the place where the gas velocity is equal to km s~^, 
which occurs at the cluster center, to the sound speed, Csc, that should occur at the cluster 
surface. Out of the cluster, pressure gradients would allow the wind to reach its terminal 
speed {v^ ^ 2csc)- 

The observed parameters of SSCs (Ho 1997; Whitmore, 2003; Turner ct al. 2003, 2004; 
Pasquali ct al. 2004; Mclo ct al. 2005; Walcher et al. 2005; Martin Hernandez ct al. 2005; 
Thompson et al. 2006; Larsen, 2006 and references therein) led us however, to re-analyze the 
original adiabatic model of Chevalier & Clegg (1985) and reahze that in the case of massive 
and compact star clusters, radiative cooling may crucially affect the internal structure and 
the hydrodynamics of the flow (Silich et al. 2004). In a more recent communication, Tenorio- 
Tagle et al. (2007) have shown that if the mechanical luminosity of a cluster exceeds a critical 
value (the threshold line), then a bimodal hydrodynamic solution for the matter reinserted 
by winds and SNe is established. In such cases, the densest central regions within the cluster 
volume undergo strong radiative cooling, what depletes the energy required to participate 
in the cluster wind. This also moves the stagnation radius, Rgt, out of the cluster center 
and promotes the accumulation and the re-processing of the enclosed matter, into further 
generations of stars (Tenorio-Tagle et al. 2005). Meanwhile, the matter injected between 
the stagnation radius and the cluster surface is still able to drive a stationary wind. This 
however, cools rapidly at a short distance from the cluster surface and becomes exposed 
to the UV radiation escaping the cluster. Wiinsch et al. (2007) has also shown that the 
location of the threshold line that separates clusters evolving in the bimodal regime from 
those with the stagnation point at the star cluster center, which are thus able to eject 
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all of the deposited mass out of the cluster, can be well approximated by simple analytic 
expressions. This is also the case for the position of the stagnation radius, the radius that 
defines the amount of matter accumulated and ejected from a cluster in the bimodal regime. 
Making use of these relationships and our semi- analytic code (Silich et al. 2004), here we 
propose another interpretation of Smith's et al. (2006) observations of M82-A1, based on 
our bimodal solution. Our model assumes that M82-A1 is embedded into a high pressure 
ISM, the thermal pressure in the core of M82 {P/k ~ 10^ cm~^ K), which is much larger 
than that found in the disks of normal galaxies {P/k ~ 10^ cm~^ K, see, for example, Slavin 
& Cox, 1993). Our final suggestion is that M82-A1 is a good example of a cluster evolving 
in the catastrophic cooling regime and that a high pressure ISM is what rapidly confines its 
wind, leading to its associated compact HII region. 

The paper is organized as follows. In section 2 we formulate our model, discuss the input 
physics and some approximations used in section 3 to develop a set of analytic equations 
which allow us to derive an approximate value of the heating efficiency required by the 
parameters of the M82-A1 cluster and its associated HII region. We confirm the analytic 
result by means of semi-analytic calculations in section 4 where we also thoroughly discuss 
the structure of the outflow driven by M82-A1. In section 5 we discuss the possible impact 
of M82-A1, if immersed into different ISM environments. In section 6 we summarize and 
discuss our results. 



2. M82-A1 embedded into a high pressure ISM 

The properties of M82-A1 and in particular its mass and expected mechanical energy 
input rate and the size and mass of its associated HII region, lead almost unavoidably to 
the conclusion that the cluster must be embedded into a high pressure region which has 
managed to confine the cluster wind. This implies a rapid hydrodynamic evolution which 
has caused the decay of the outer shock, inhibiting the formation of an outer shell of swept 
up matter and the built up of a superbubble (as otherwise described in Weaver et al. 1977). 
Instead, a standing reverse shock sitting near the cluster surface, decelerates and thermahzes 
the outflow increasing its thermal pressure. The shocked gas then slows down while it cools 
and is displaced away from the cluster, as fresh wind matter enters the shock. The shocked 
gas thus recombines at a small distance behind the reverse shock and becomes exposed to 
the UV radiation from the cluster, forming a stationary HII region. Note that upon cooling 
the shocked wind gas ought to undergo a rapid condensation to counterbalance the lost of 
temperature (from the post-shock temperature to the assumed photoionization temperature 
Thii ~10^ K) and be able to keep a pressure similar to that of the surrounding ISM (Pism)- 
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In this model the HII region associated to the central cluster results from the photoionization 
of three different regions. The central zone within Rgt where the injected matter cools in a 
catastrophic manner. A second photoionized zone is generated within the free wind region 
as this gas cools suddenly to Thu ~ 10*^K and finally, the shocked wind gas that has more 
recently being able to cool down from the high temperatures acquired after crossing the 
reverse shock. Other wind material that experienced earlier the same evolution is steadily 
pushed out of the standing ionized outer shell by the newly incoming material, what allows 
for its recombination, its further cooling by radiation and to be further condensed so that 
its pressure remains equal to Pism- 

Major restrictions to the model arise from the pressure balance between Pjsm and the 
wind ram pressure. This defines the position of the standing reverse shock, the size of the 
coohng distance that the shocked gas ought to travel before becoming exposed to the cluster 
UV radiation and which together with the available number of photons, defines also the 
size of the outer HII region, which must agree with the observations. Here we show that to 
set Rhii in agreement with the observations, the location of Rah and the size of its cooling 
distance imply that M82-A1 must be experiencing a bimodal hydrodynamic solution. In 
the bimodal regime, strong radiative cooling largely diminishes the cluster thermal pressure 
within the densest central zones. This leads to the accumulation of the matter reinserted 
within the cooling volume, reducing the amount of mass that a cluster returns to the ISM. 
This also results into a slower acceleration of the matter injected in the outer segments of 
the cluster (r > i?st), what leads to a reduced terminal speed of the resultant wind (see 
Tenorio-Tagle et al. 2007; Wiinsch et al. 2007). As shown below, these considerations allow 
us to predict which is the heating efficiency, or the fraction of the total energy input rate that 
is used to produce the wind of M82-A1 and which is the structure acquired by the outfiow. 

To derive an analytic expression for the position of the stagnation point, Wiinsch et al. 
(2007) used the mass conservation law together with the fact that the density at the star 
cluster surface scales almost linearly with that at the stagnation point. They found that for 
star clusters with a given radii Rsc-, the location of the stagnation radius depends on the 
excess star cluster mechanical luminosity above the critical value, I/cr«t- 
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where < < 1 is the heating efficiency or fraction of the deposited kinetic energy that 
after full thermalization is not immediately radiated away, and thus is instead evenly spread 
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Fig. 1. — Schematic representation of the M82-A1 HII region. The ionizing cluster is con- 
tained within the radius Rsc- Rst is the stagnation radius that separates the inner cluster 
region, where catastrophic cooling sets in, from the outer zone {Rst < r < Rsc) where the 
thermalized ejecta remains hot and is thus able to compose a stationary wind. The star 
cluster wind however cools rapidly and becomes exposed to the UV radiation from the clus- 
ter forming a high speed (Voo) HII region component. Eventually this gas is stopped and 
re-heated by the reverse shock which sits at a fixed radius, Rgh, defined by the condition 
that the wind ram pressure equals that of the ISM. The wind processed at the reverse shock 
then cools down at a distance -Ra to become once again target of the cluster UV radiation, 
causing in this way the outer, photoionized, thin standing shell with radius Rhii (all scales 
are greatly distorted). The model assumes an even pressure in all regions: from the ISM 
to the location of the reverse shock where the wind ram pressure equals that of the ISM. 
Regions filled with photoionized material are shaded. 
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within the star cluster volume, causing the overpressure that drives the wind (Stevens & 
Hartwell 2003, Melioli & de Gouveia Dal Pino 2004). Here a — 0.28 is a fiducial coefficient 
(see Wiinsch et al. 2007), jii = 14m ///II is the mean mass per ion, Vaoo = C^Lsc/ MscY^'^ 
is the adiabatic wind terminal speed, Msc is the mass deposition rate, Cst and Agt are the 
speed of sound and the value of the cooling function at the stagnation point, respectively. 
Note also that the density and the temperature at the stagnation radius are related by the 
equation (see Silich et al. 2004; Wiinsch et al. 2007): 



rist 



yvlj'^ 4/(7 - 1) 



nl/2 



(3) 



where qm = (3M5C')/(4/t7?|(^) is the mass deposition rate per unit volume. 



Figure 1 presents a schematic representation of the model. M82-A1 is contained within 
the inner solid line. The dashed line interior to this marks the location of the stagnation 
radius. Interior to Rgt, cooling dominates over the heating provided by stellar winds and 
supernovae. The injected matter there cools down rapidly and is unavoidably reprocessed 
into new generations of stars supporting a low level of star formation as the parent cluster 
evolves (see Tenorio-Tagle et al. 2005). 

The outer cluster zone, Rgt < r < Rsci is filled with hot thermalized ejecta unable to 
radiate away all of its thermal energy within a dynamical time scale, and thus capable of 
moving out of the cluster composing a stationary wind. The amount of mass reinserted by 
the cluster wind is (Wiinsch et al. 2007): 
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Outside the cluster there is a free wind region {Rsc < r < Rsh) bound by the standing 
reverse shock at a radius Rgh- The free wind is also able to cool rapidly, loosing its thermal 
pressure while becoming easy target of the cluster UV radiation field and thus forming 
another component of the HII region associated to the cluster. At the reverse shock the 
ejected material, which streams with its terminal velocity V^, is decelerated and re-heated. 
The matter behind Rgh would remain hot until it reaches the cooling radius R\ 



Ra — Rsh + L\, 
where the coohng length La is (Franco, 1992): 
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where He = 14mH/23 is the mean mass per particle in a fully ionized plasma that contains 1 
helium atom per every 10 atoms of hydrogen, k is the Boltzmann's constant, T, = 3fieV^/16k 
and Ug = (7 + l)/(7 — l)p^/fii are the post-shock temperature and the post-shock ion 
density, respectively, A(Ts) is the cooling function and fx is a fiducial coefficient that takes 
into consideration the enhance cooling promoted by the continuous growth of density in the 
post-shock region. Good agreement between formula ([6]) and numerical calculations was 
found for fx = 0.3 (see below). 

Outside i?A, the wind gas is photoionized by the UV radiation produced by the star 
cluster, forming a standing photoionized shell, with a pressure identical to that of the ISM. 

Note that the location of the various disturbances depends only on the assumed value 
of PisM and on the terminal speed of the cluster wind. In this way, the density of the wind 
at the shock radius, p^, is 

Pw = Pram I V^, (7) 

where Pram, the wind ram pressure, is equal to Pism- In the strongly radiative case the 
terminal speed of the wind, V^, falls below the adiabatic value. Nevertheless, one can use 
the adiabatic relation between V^o and c^t, the sound speed at the stagnation point (see 
Canto et al. 2000), as a good approximation that allows one to determine the radiative wind 
terminal speed if the sound speed at the stagnation radius is known: 



where the radiative energy losses across the cluster were assumed to be identical to those at 
the stagnation point. The effective energy deposition within the cluster is then: 

2 

qeff = me - n^Ast = -Trf- Qe, (9) 

1 — ^ ^Aoo 

where Qe is the average energy deposition rate per unit volume within the cluster. 



3. The analytic model 

The structure of the outflow can be derived analytically from a set of equations that 
consider: The available number of ionizing photons, photoionization balance, pressure con- 
flnement and the divergency of the wind stream: 
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where f3 = 2.59 x 10^^^ cm^'^ s^^ is the recombination coefficient to all but the ground level, 
N^'^ is the total number of UV photons produced by the cluster and A^°"* is the total number 
of photons escaping the cluster and the free wind region, and thus available to impact on the 
shocked wind gas once this has cool by radiation and has acquired by condensation a density 
uhii- The stagnation radius, Rgt, and the amount of matter that a cluster returns to the ISM, 
Mout, are defined by equations ([1]), ([2]) and (jl]). /t is the fraction of the star cluster ionizing 
radiation that reaches the cool shocked gas behind the standing reverse shock. Combining 
equations (fTTI) . and ([5]), one can derive the expression for the reverse shock radius, Rsh- 
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On the other hand, the position of the reverse shock is defined by the mass conservation (fT3l) 
and pressure balance (fT2!) equations and by the amount of mass that the cluster returns to 
the ISM (equation m): 
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where \4o and Vaoo are the radiative and adiabatic wind terminal speeds, respectively (see 
Wiinsch et al. 2007). 

Substituting equation (fT5|) into equation (fT^ we obtain a nonlinear algebraic equation 
which defines the heating efficiency, t], required to match the observed parameters of the HII 
region {Rhii, N"^"^) associated to M82-A1, if one accounts for the star cluster parameters 
{Rsc, Lsc) and the inferred pressure of the ISM. 
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If the adiabatic wind terminal speed, Vaoo, is known, the only free parameter in equation 
f|T6|) is < < 1, the fraction of UV photons able to reach the outer HII region. One can 
then solve equation ( fT6|) by iterations for a given value of ft- 
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The solution of equation f|T6|) shows that the heating efficiency which is required in order 
to fit the observed structure of the M82-A1 HII region is small, 1] ^ 4.65%, and does not 
depend significantly on the adopted value of ft- Figure 2 displays the structure of the outflow 
produced by M82-A1 as well as the dimensions of its associated outer HII region, calculated 
under the assumption that V^oo = 1000 km s~^. The low value of the heating efficiency 
results into a large stagnation radius {Rgt ~ 2.94pc, the lower solid line in Figure 2) which 
implies also that a large fraction of the matter supplied by stellar winds and supernovae is 
to remain bound to the cluster and thus only the small amount deposited between Rgt and 
Rsc would conform the cluster wind while expanding approximately with 200 km s~^. The 
outflow is strongly decelerated and re-heated at the reverse shock (4.24pc < Rgh < 4.35pc, 
the second solid line from the bottom in Figure 2) and then cools down rapidly to form 
an outer shell of photoionized matter depleting the photons that escape the cluster (A^°"*). 
Note that the thickness of this outer shell is very small, and becomes even smaller as one 
assumes a smaller number of photons escaping the SSC and the cluster wind. 



4. The semi-analytic model 

To corroborate our results, here we use our semi-analytic code (see Silich et al. 2004) that 
accurately calculates the divergent spherically-symmetric outflow beyond Rst, accounting for 
strong radiative cooling. 

We use ffist the trial heating efficiency r] to calculate the position of the stagnation 
point as it was suggested in Tenorio-Tagle et al. (2007) and then integrate the equations of 
mass, momentum and energy conservation outwards, knowing that at some distance from 
the cluster the standing reverse shock sets in. We stop the integration at the reverse shock 
and use then the Rankine-Hugoniot conditions to calculate the thermal pressure and the 
velocity of the plasma behind the shock front: 



V 7 + 1 

7-1 2 1 
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U2 = U^( + —^T72 ' (18) 



where and u^, and P2 and U2 are the thermal pressure and the velocity of the outflow 
ahead and behind the reverse shock, respectively. Mi = Uw/cw is the Mach number ahead 
of the shock front, and 7 = 5/3 is the ratio of specific heats. P2 and U2 together with the 
mass conservation law, Mout = 47rp^-u^i?^^, define the initial conditions behind the reverse 
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Fig. 2. — The analytic model of M82-A1 and its associated HII region. The various hnes 
mark the location of different hydrodynamical disturbances as a function of /<. Prom bottom 
to top these represent the location of the stagnation radius (Rst), the cluster radius (Rsc): 
the reverse shock location (Rsh), the cooling radius, R\ and the outer radius of the standing 
HII region {Rhii)- 
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Fig. 3. — The semi-analytic model of M82-A1 and its associated HIT region. Panels a - d 
present the distributions of the velocity, thermal (solid line) and ram (dotted line) pressures, 
density and temperature, respectively. The size of the standing outer HII region is indicated 
in panels c and d by the two vertical dotted lines. The calculations assumed an ft — 0.5, 
Z — Zq and an adiabatic wind terminal speed, Vaoo — 1000 km s~^. 
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shock for the set of main equations outside of the cluster (see Sihch et al. 2004): 

du^ _ 1 (7 - l)rQ + 2'-fUu,Pw 
dr Pw r{ul, - cl) 

dPyj _ Msc duu, 
dr Anr'^ dr ' 



(19) 
(20) 



Msc 

= -A 2- (21) 

The run of the hydrodynamical variables in the outer part of the flow was obtained by 
integrating equations (fT9l) - (l2Tl) from the reverse shock radius outwards. 

The thermal pressure, P^, and the wind expansion velocity, Uw, ahead of the shock 
front depend on the star cluster parameters and on the radius of the reverse shock, Rgh. 
Therefore the set of initial conditions f|T7|l - f|T8l) contains two model parameters: the value 
of the heating efficiency, 77, and the position of the reverse shock, Rgh- We iterate rj and Rsh 
until the conditions 



Ra — Rhii 



1 - 



3N 



out 



1/3 



(22) 



Phi I = PisM, (23) 

are fulfilled. Here Rhii = 4.5pc is the observed radius of the M82-A1 HII region, Prii = 
kriHiiTHii and Thii = lO^K are the thermal pressure and the temperature of the ionized 
gas in the standing outer shell, respectively. 

Figure 3 presents the results of the calculation assuming ft = 0.5, e.g. when only 
half of the ionizing photons reach the outer standing shell. The expansion velocity of the 
outfiow grows from zero km s~^ at the stagnation radius to rapidly reach its terminal value, 
Voo ~ 180 km s~^, and then drops when the free wind reaches the reverse shock at r ^ 4.4pc 
(panel a in Figure 3). The gas decelerates rapidly at the reverse shock and then as it cools 
down it is further condensed. Thus in the case of the pressure confined wind, the reverse 
shock separates gas that flows supersonically in the free wind region from gas moving away 
with a subsonic velocity. The thermal pressure (see panel b) drops initially when the gas is 
accelerated to reach its terminal speed, and remains at its lowest value as a consequence of 
radiative cooling. Pt^ is restored after crossing the reverse shock, where the ram pressure 
in the free wind reaches a balance with the thermal pressure in the surrounding interstellar 
medium, and finally after photoionization of the dense outer layer. The density decreases as 

in the free wind region. It is then compressed at the shock, reaches, upon condensation 
induced by radiative cooling, the HII region value and finally the maximum value in the 
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neutral layer outside of the outer photoionized shell (Panel c in Figure 3). Figure 3 also 
shows that in the bimodal regime the temperature of the outflowing matter drops rapidly 
outside of the cluster. The free wind is then re-heated at the reverse shock and cools down 
again to form the outer shell whose inner skin is photoionized by the Lyman continuum 
escaping from the M82-A1 cluster. Note that in our calculations the temperature of the free 
wind (between the star cluster surface and the reverse shock) drops to lO^K and thus a broad 
{uw ~ 180km s~^), low intensity emission line component is expected from this region. 

Certainly, the analytic model is unable to show the detailed structure of the outflow. 
Nevertheless, the analytic heating efficiency, stagnation radius, reverse shock and cooling 
radii are in reasonable agreement with the results obtained from the semi-analytic model: 
r] = 6.8%, Rst = 2.86pc, Rsh = 4.41pc and R\ = 4.45pc, if the fraction of ionizing photons 
reaching the outer HII region is ft = 0.5. 



5. M82-A1 in different environments 

The theory developed in the previous sections is based on the assumption that the 
HII region detected around the M82-A1 is a standing, photoionized shell of shocked wind 
matter, confined by the high pressure of the ISM. A major implication of the results above 
is that clusters with properties (size, mass and age) similar to those of M82-A1 would have 
a low heating efficiency. Here we assume such a cluster to be embedded into different ISM 
environments and workout the properties of the resultant pressure confined remnant. 

The different pressure in the ambient ISM does not affect the distribution of matter 
within the cluster. It would not affect the position of the stagnation radius, nor the density, 
temperature and velocity of the ejected plasma. A different ISM pressure would modify only 
the outer structure of the outflow shifting the position of the reverse shock. This is because 
of the pressure equilibrium condition: Pram = Pw{Rsh)V^ = Prii = Pism- Using equations 
dl]) and f[T^ one can then obtain the position of the reverse shock for M82-A1 sitting in 
different interstellar environments: 

yAirPisM J yA'KPisM ) \Lsc J 

One can then flnd the outer radius of the HII region from equations ([5]), ([6]) and (ITT]) , knowing 
that the density in the ionized outer shell is a linear function of the pressure in the ambient 
interstellar medium: 

riHii = PisM/kTHii, (25) 
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where the temperature of the photoionized shell, Thii^ is set through photoionization and 
therefore it does not depend on Pism- Equation fl2^ indicates also that there exists the 
critical interstellar pressure, Pent'- 



If the pressure in the ambient ISM exceeds this critical value, Pism > Pcru, the cluster would 
not have the sufficient power to drive a cluster wind. In this case all matter deposited by 
stellar winds and SNe would remains buried within the cluster. 

Figure 4 shows the structure of the pressure confined remnants that M82-A1 would 
produce in different interstellar environments. Our results show that the outer ionized shell 
is very thin. Also that the radius of the ionized shell shrinks towards the star cluster if one 
considers a pressure in the interstellar medium that approaches the critical value (equation 
l26i) and grows as the considered pressure in the interstellar medium drops. The analytic 
results (dashed and dotted lines in Figure 4) show reasonable agreement with those obtained 
by means of semi-analytic calculations (cross symbols in Figure 4). 



The results from our analytic formulation to a pressure confined wind in agreement with 
those obtained with our semi-analytical code, imply that M82-A1 is a massive and compact 
cluster with a low heating efficiency. This implies a completely different result from what 
one would expect from an adiabatic model, as a low heating efficiency leads to a bimodal 
hydrodynamic solution and with it to a low mass deposition rate into the ISM with a much 
reduced outflow velocity. 

Furthermore, to match the observed parameters of M82-A1 and its associated HII region, 
our results lead also to a high pressure environment able to conflne the cluster wind by setting 
a reverse shock close to the star cluster surface. In this way the outflow is thermalized, what 
leads to temperatures near the top of the interstellar cooling curve and thus to a rapid 
cooling of the strongly decelerated outflow. The wind becomes then target of the cluster 
UV radiation, composing a narrow standing outer shell of photoionized gas with the shocked 
wind matter that continuously traverses the reverse shock. 

Our calculations lead to three HII region components associated to M82-A1: A central 
component where the deposited matter cools catastrophically and does not participate in the 
wind. A section of the free wind region which, upon expansion, is also able to cool rapidly. 
And the outer stationary shell of shocked wind matter, which deflnes the observed size of 




1/2 



(26) 
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Fig. 4. — The structure of pressure confined winds into different ISM environments. We 
have assumed the star cluster parameters derived by Smith et al. (2006) for M82-A1 and a 
heating efficiency, rj = 4.6%, as derived in section 3. The sohd fine marks the star cluster 
radius which is equal to the half-light radius of M82-A1. The dashed and dotted lines mark 
the position of the reverse shock, Rgh, and of the outer Hll region, Rhii, respectively, as a 
function of Pjsm normalized to that found by Smith et al. (2006) in the M82-A1 Hll region. 
The cross symbols represent the results of the semi-analytic calculations and also mark the 
positions of the shock front and the outer edge of the HII region (lower and upper symbols, 
respectively) . The fraction of the star cluster ionizing radiation able to photoinize the outer 
shell is ft — 0.5 in both, the analytic and semi-analytic calculations. 
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the associated HII region. Note that, given the densities in the latter component, this ought 
to be the one producing the most intense emission hues, while the second component is to 
produce a broad ( 2 Vqo) low intensity emission line component, similar to that detected 
in NGC 4214-1 (sec Ho & Fihppenko 1996), able to enhance, as in the case of M82-A1, the 
width of the observed lines. 

Our two methods of solution have also shown good agreement when considering the 
wind of M82-A1 being pressure confined by an ISM with different pressures. 

Note that in the case of very young SSCs (a few Myrs old) their mass is usually derived 
from the intensity of their optical emission lines, emanating from their associated compact 
HII regions and thus through a measure of the available UV photon flux, an assumed IMF 
and stellar synthesis models (Melo et al. 2005, Smith ct al. 2006). In the light of the models 
here presented, this implies both that young SSCs are embedded into a high pressure ISM 
and also that their heating efficiency ought to be rather small for their winds to be pressure 
confined in the immediate neighborhood of the clusters, and thus present an associated low 
mass and compact HII region. 

Note also that the two different populations of massive and compact clusters (Gyr old 
and very young clusters) present a similar range of masses as well as a similar size distribution. 
This can be explain if most of the mass reinserted by stellar winds and supernova explosions 
during the early evolutionary stages (which may amount to 30 per cent of the Msc released 
during the first 50 Myr of evolution) is not returned to the ISM but rather re-processed in 
situ into further stellar generations. This would allow SSCs to keep most of their mass (and 
size) as they age and are able to avoid their dispersal. This also implies that a large fraction 
of the low mass stars presently observed in old compact and massive stellar clusters might 
have formed from the matter injected by massive stars. In this respect our bimodal model 
provides the hydrodynamical ground to a self-enrichment scenarios. In these, for example, 
the abundance anomalies observed in globular clusters may result from the enrichment of 
the proto-stellar gas during the earher stages of their evolution (see, for example, Prantzos & 
Charbonnell, 2006; Decressing et al. 2007 for a comprehensive discussion of different aspects 
of this problem) . Our solution also suggests that the metallicity of the compact HII regions 
associated with the intermediate age (5 Myr - 10 Myr) massive clusters should be super-solar 
and variable with time. 

We thank our anonymous referee for valuable comments and suggestions. This study 
has been supported by CONACYT - Mexico, research grant 47534-F and AYA2004-08260- 
C03-01 from the Spanish Consejo Superior de Investigaciones Cientificas. 
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